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Colloidal particles or nanoparticles, with equal affinity for two flu- 
ids, are known to adsorb irreversibly to the fluid-fluid interface. We 
present large-scale computer simulations of the demixing of a bi- 
nary solvent containing such particles. The newly formed interface 
sequesters the colloidal particles; as the interface coarsens, the par- 
ticles are forced into close contact by interfacial tension. Coarsening 
is dramatically curtailed, and the jammed colloidal layer seemingly 
enters a glassy state, creating a multiply connected, solid-like film 
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in three dimensions. The resulting gel contains percolating domains 
of both fluids, with possible uses as, for example, a microreaction 
medium. 

The search for new materials with mesoscale or nanoscale structure is a major theme 
of current physical science. Routes which exploit spontaneous self assembly in thermal 
equilibrium are certainly important PIE], but nonequilibrium processes offer more control 
- since assembly is then governed not just by thermodynamic conditions but by the entire 
process history (e.g. OElll]). Moreover, the resulting materials may become trapped 
in deeply metastable states such as colloidal clusters, glasses, and gels [HI El El EI, re- 
maining more robust than an equilibrium phase to subsequent changes in thermodynamic 
conditions. This is a key consideration in determining the shelf-life and flow behavior of 
everyday products such as paint, vaccines, and yoghurt. 

We use computer simulations to explore a kinetic pathway that leads to the creation 
of amorphous soft-solid materials. In such a material, which we call a 'Bicontinuous 
interfacially jammed emulsion gel' (or 'Bijel', for short) [S], a pair of interpenetrating, 
bicontinuous fluid domains are frozen into a permanent arrangement by a densely jammed 
monolayer of colloidal particles at the fluid-fluid interface. Such materials may have some 
remarkable properties, stemming directly from the nonequilibrium, arrested nature of the 
monolayer: Bijels should be highly tunable in elasticity and pore size through the volume 
fraction and radius a of the solid particles. (The radius can be varied from microns 
to nanometres without altering the physics of structure formation by the route reported 
here.) One possible application, explored below, is as a cross-flow microreaction medium 
in which two immiscible fluids are continuously brought into intimate contact by pumping 
them in opposite directions through a static Bijel. 

To achieve maximal stability of a particle-laden interface, the colloidal particles should 
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be chosen with nearly equal affinity for the two liquids involved jH]. This creates simi- 
lar values for the two fluid-solid interfacial tensions, and thus a fluid-fluid-solid contact 
angle close to 90 degrees (known as 'neutral wetting'). A spherical particle is then in 
stable equilibrium with its equator positioned at the fluid-fluid interface. In practice, this 
equilibrium is so stable that detachment of such a particle cannot be achieved by thermal 
motion alone jH] . For neutral wetting, the fluid-solid interfaces have the same total energy 
regardless of particle position, but the fluid-fluid interfacial area is reduced, by a disc of 
radius a, when the particle lies midway across the interface. The detachment energy e 
is the interfacial energy of this disc, e = crvra^, with a the fluid-fluid interfacial tension. 
Hence e/ksT = (a/ao)^ where = kBT/na. For T = 300 K and typical a of order 0.01 
Nm~^ or larger, Oq is 0.4 nm or less. Thus e/ksT > 10 even for a particle of 1 nm radius, 
and thermally activated detachment can be safely neglected for, say, a > 3 nm. 

We consider near-neutral wetting particles suspended in a binary solvent under con- 
ditions where the fluids are fully miscible (generally at high temperature) and of roughly 
equal volume fraction. In the absence of strong attractions between them, the particles 
will diffuse freely. If the temperature is quenched deep into a two phase region, the sol- 
vents will demix by spinodal decomposition ^U]. A sharp interface between the two fluids 
develops, and coarsens. During the coarsening, which is driven by the tendency of the 
interface to reduce its area, the characteristic lengthscale L{t) initially increases with time 
in a well-understood manner ^I], causing bumps on the interface to flatten, and causing 
necks between neighboring domains of the same fluid to pinch offi 

We have studied what happens after this initial separation by simulations ^21 using the 
lattice Boltzmann (LB) method [ini CH CHI IIEI El • We find that, as coarsening proceeds, 
the interface sweeps through the fluid phases, efficiently collecting and sequestering the 
colloidal particles. Initially, the attached particles have little effect on the interfacial 
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motion, but as more are collected and the interfacial area shrinks, they soon approach a 
densely packed monolayer. At this point, the fluid must either (i) stop coarsening at some 
lengthscale L{t) = L* or (ii) thereafter expel particles steadily from the ever- shrinking 
interface. In our simulations, we see a drastic curtailment of the coarsening and very little 
particle expulsion. This suggests that the free energy landscape of the dense colloidal film 
is such as to trap particles on the interface in a metastable, amorphous state fHl- 

Fig.l shows the structure as time evolves; for visibility, only a small part (in cross- 
section, upper panels, or crop, lower panels) of the full system is shown. Ref.jSOj shows 
an animation of the observed sequence of events within the cropped region of Fig. 1; this 
clearly shows the particle sweep-up and the dramatic slow-down of coarsening. Ref. |21j 
shows the entire simulation domain for another run, in which bidisperse particles are 
used; this prevents development of any local crystalline order in the interfacial colloid 
layer (arguably visible in Fig.l). Ref. [22] shows a similar subregion for this case as in 
Ref. [201 • The results are qualitatively the same: see Fig.2 for a snapshot at late times. 
Fig.3 shows the time dependence of the domain size L{t) in each case. 

The parameter values chosen for these simulations map onto particles of radius 
a = 5 nm in a symmetric pair of fluids each having viscosity rj = 10^^ Pa s and mass 
density p = 10^ kg m~^ , with cr = 6 x 10"^ Nm^^ at T = 300 K; such values are 
typical of a short-chain hydrocarbon/water or alcohol/water mixture. Our particles have 
purely repulsive interactions, with range extending somewhat beyond their hard-sphere 
(hydrodynamic) radii ^2!, so that particles remain visibly separated even in a dense 
monolayer. The parameter mapping is made by matching dimensionless control groups 
e/kT and apa/yf [23 E]- Brownian motion of the colloidal particles is included ^2], but 
has rather little effect during the time regime we can reach by simulation (see below), and 
would have even less effect with larger particles. We have also checked the role of short- 
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range, thermally reversible bonding among colloids [7j, but this too has little effect. These 
observations are attributable to the strong separation between Brownian and interfacial 
energy scales referred to above. 

A number of numerical compromises were made to keep the simulations tractable ^2] • 
First, the Reynolds number Re = [dL / dt) pa / rj is much larger than in the real system, 
though we still have Re <C 1 [SHI- Second, the scale separation between the particle radius 
a and the fluid-fluid interfacial thickness ^ is only modest (a factor two or three), with 
the lattice spacing, in turn, not much less than ^. This gives imperfect discretization 
of the fiuid-fiuid interface around particles and may under-represent the energy barrier 
to short-scale rearrangements ^2]- Finally, for the physical parameters given above, the 
effective run-time of our largest simulations is about 300 ns. (For larger particles, say 
a = 3/im, the equivalent run time would be around 5 ms.) 

Although these simulations dramatically confirm our proposed kinetic pathway for 
creating a fiuid-bicontinuous state with a particle-laden interface, they cannot tell us 
whether this state is a fully arrested gel on laboratory time scales. Unlike our simulations, 
the latter vastly exceed the time scale r = QTirja^ /ksT characterising Brownian motion 
of colloids. Nonetheless, the observed behavior, particularly for bidisperse particles, is 
consistent with that of a long-lived, metastable, arrested state. In common with other 
such states (including colloidal glasses), Bijels might show slow aging behavior as the 
interface approaches saturation. This may explain the residual late-time dynamics visible 
in the L{t) curves of Fig.3. Alongside aging, the slow residual dynamics could be due to 
the incomplete separation of length scales in LB noted above or, in the monodisperse case, 
due to a tendency for the interfacial layer slowly to acquire local crystalline order. (Such 
ordering would not preclude, and might even enhance, eventual structural arrest.) We 
have also assessed the particle mobility in the interfacial film by measuring the distribution 
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of individual particle displacements at late times. We found this to be dominated by the 
residual relaxation dynamics of the structure rather than by diffusion within the film. 

The near-complete cessation of fluid motion on time scales of order r suggests that 
further coarsening, which requires expulsion of particles from the interface, cannot take 
place without colloidal Brownian motion. If, as argued earlier, such motion is ineffective at 
detaching particles from a static interface, coarsening must cease altogether. In principle, 
lateral diffusion within a film of fixed area might continue. However, the surface pressure 
of such a film is estimated ^Hj as 11 ~ ksT/la'^^ipm — V^)] where ip is the areal fraction 
and ipm that of a maximally close-packed configuration. For this pressure to balance the 
interfacial tension requires ipm — ip ^ {ao/ay. For large enough particles, this ensures 
that ip exceeds any threshold value, so long as this value is less than -i/^m, for the onset of 
an arrested state within the film. 

Further evidence for interfacial arrest was gained by additional, higher-resolution sim- 
ulations that examine the dynamics of two specific structural motifs characteristic of the 
bicontinuous structure. One of these is a long cylinder, representing a fluid neck. With- 
out particles the Rayleigh-Plateau instability j23] would cause the cylinder to break into 
droplets. We show in Fig. SI a dense bidisperse colloidal packing on a cylindrical interface. 
When perturbed, this shows no sign of either ordering or breakage, and the initial undula- 
tion visibly decays, rather than grows. The structure then arrests, and remains unchanged 
for the duration of the simulation, which is ten times longer than the time-to-rupture of an 
unprotected cylinder. Our second structural motif is a periodically rippled surface, whose 
bumps are broadly characteristic of a non-necklike section of the bicontinuous interface 
between fluids. Without particles, the ripple would rapidly be pulled flat by interfacial 
tension. Fig.S2 shows how this process is interrupted by interfacial jamming. After an 
initial transient decay, the film jams, and bumps on it persist at least thirty times longer 
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than the decay time without particles. There is neghgible macroscopic motion during the 
latter half of this simulation. 

These results show that, at sufficient interfacial coverage both necks and bumps 
can arrest by jamming of the adsorbed colloidal layer into a densely packed, and therefore 
presumably glass-like, state. Since these two structural elements are, in combination, 
the driving features of bicontinuous coarsening [21], their arrest would be enough to 
prevent coarsening. Hence these studies provide very strong supporting evidence for an 
eventual arrest of our bicontinuous structure, caused by a jamming transition [2^] within 
the colloidal monolayer. This transition is induced by the (capillary) stresses arising from 
the fluid-fluid tension. 

Once the interfacial film has fully arrested, since it percolates in three dimensions, the 
entire material will acquire solid elasticity at scales beyond L*. The static modulus G of 
the resulting gel should scale with the interfacial energy density cr/L; so long as nearly 
all particles end up on the interface, L* ^ a/(j) with (p the particle volume fraction. For 
a = 0.01 Nm-\ 0.01 < < 0.1 and 5nm < a < 5/im we estimate 20 < G < 2 x 10^ 
Pa. This is a very wide 'tuning' range for material design. Under nonlinear stress the 
interfacial area will dilate significantly: only a modest dilation (say 10%) may suffice to 
cause melting of the particle layer and drastic fiuidisation. This might instigate both fiow 
and coarsening above some yield stress Y ^ O.IG. If the stress falls back below Y, we 
expect resolidification, possibly with remanent anisotropy (hysteresis). The nonlinear fiow 
behavior of these gels could thus show a remarkable strain-melting, possibly reminsicent 
of a colloidal glass, but with a much higher stress scale set by interfacial, not Brownian, 
forces. The estimates above for the material properties of the gel stem from the jamming 
of colloids by the interfacial forces and apply even for purely repulsive particles. Any 
additional bonding attraction, if of sufficient strength, might enhance the rigidity of the 
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interfacial layer, but may also cause colloidal aggregation within the bulk phase(s) prior 
to monolayer formation jSZj. Fusing the colloids after gel formation (e.g. by irradiation) 
would completely stabilize the structure and drastically alter the flow behavior. 

Our study differs from previous work in which colloidal particles were used to stabi- 
lize spherical jHl EH] and aspherical j2Zj emulsion droplets of one liquid in another. Under 
compression, such emulsions can form robust gel phases with interesting mechanical prop- 
erties j2n], but fluid bicontinuity is not among them. The previously preferred route to 
particle-stabilised emulsions involves agitation of immiscible fluids and does not appear to 
favour bicontinuity P . Other related work jHl EDI 1^ involves particles with a strong 
preference for one of the two liquids, creating a particle network within the chosen liquid 
rather than at the interface. 

We expect Bijels to have several interesting physical properties beyond those discussed 
above. First, the fluid-bicontinuous state should remain equally insoluble on exposure to 
either of its solvents. This contrasts with particle-stabilized emulsion gels formed by 
compression [201, in which an excess of the continuous phase could cause droplets to 
separate, losing macroscopic rigidity. In Bijels this will not happen, since neither of 
the two interpenetrating fluids can alter its volume without also increasing the total 
interfacial area. The Bijel can thus metastably support simultaneous coexistence with 
bulk phases of both fluids. This is reminiscent of an equilibrium property of middle-phase 
microemulsions, which are bicontinuous states with an interface stabilized by amphiphilic 
molecules jHB]- In contrast to Bijels, such microemulsions are not gel phases, but inviscid 
fluids, result of their high interfacial mobility. 

Second, fluid bicontinuity imparts high permeability of the gel to either of its compo- 
nent solvents, and any reagents dissolved in them. Accordingly Bijels may have potential 
as media for continuous-process microreactions [S^ ES]. Specifically, a static gel could 
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support a steady permeation flow of botli fluids simultaneously in opposite directions. 
This would bring two molecular reagents, soluble only in mutually immiscible fluids, into 
intimate contact at the fluid-fluid interface in the interstitial regions between the col- 
loids. Reaction products soluble in either phase would be swept out continuously. To test 
this concept, Fig.4 shows a simulation in which the two fluids are moving through the 
structure in opposite directions. On the timescale of our simulation, the gel has easily 
enough mechanical integrity to sustain this cross-flow without breaking up. Within the 
mapping onto physical parameters made previously, the chosen cross-flow fluid velocity 
V = 0.01cr/?7 is of order 10 cm s~^: this is an extremely large value, given the pore scale 
L* of only 70 nm. Local shear rates are of order 10^ s~^. 

In summary, we have presented simulation data showing formation of a self-assembled 
bicontinuous structure with interfacially sequestered particles. This followed a kinetic 
pathway involving a colloidal suspension in a binary solvent, initially miscible, that un- 
dergoes a temperature quench. Our simulations show a drastic curtailment of coarsen- 
ing, consistent with eventual structural arrest: a scenario further supported by higher- 
resolution studies of appropriate structural motifs (bumps and necks). This suggests a 
route to the creation of new class of gels, Bijels, with potentially remarkable physical 
properties. 
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Fig.l. Time evolution of monodisperse neutrally wetting colloidal particles at 
volume fraction = 0.2 in a binary solvent following a quench. In the upper panels, a 
128 X 128 section through the full simulation (performed on 128^ lattice) is shown. The 
two fluids are colored yellow and blue. Monodisperse colloidal particles (green) are shown 
only if overlapping the plane of the section, with parts lying behind this plane occluded (so 
that particles whose midpoint is behind the plane appear reduced in size). Sequestration 
of the particles, followed by near-arrest of the bicontinuous structure, is seen. Frames 
are at 5 000, 50 000, 500 000 LB timesteps; for parameter details see jT^j. In the lower 
panels we show a three dimensional view of a 64 x 64 x 24 piece cropped from the same 
simulation with the same time-frames. Both fluids are now shown transparent with the 
two sides of the fluid-fluid interface painted yellow and blue. Sequestration, arrest and 
fluid bicontinuity of the resulting structure are all apparent. 
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Fig. 2. Structure at 200 000 timesteps in a simulation of bidisperse neutrally 
wetting particles in a binary solvent following a quench, in section and 3D view. Color 
coding as in Fig.l but with larger particles shown in red; size ratio 1.2:1. On the left 
is a 128 X 128 section; the top left quadrant of this square is coincident with the front 
face of the 64 x 64 x 24 cuboid used for the 3D view (at right). There is no sign of 
local crystalline order in the particle film (visible in Fig.l, lower panels) although some 
tendency toward local segregation by particle size is seen. For parameter details see [T^ . 
The same structure, evolved to 590 000 timesteps, is shown in Fig. 4 below. 
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Fig. 3. Time evolution of the structural length scale L{t) for monodisperse 
(lower curve) and bidisperse (upper) particle runs as described in [T2]. L is measured in 
units of the LB lattice and t in LB timesteps. Without particles, coarsening would proceed 
with the slope indicated by the dotted line (see [H]). Near-arrest is visible, suggesting a 
finite asymptotic domain size L*, particularly in the bidisperse case. This L* is less than 
1/4 of the simulation box size and not limited by finite-size effects jTTj. For the parameter 
values chosen (corresponding to 5nm particles) the data shown run approximately from 
6 ns to 300 ns in real time jT2] , with L* ^ 70 nm. (At times less than 6 ns, the fiuids are 
demixing diffusively so that sharp interfaces have not yet formed.) 
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Fig. 4. A 64^ X 24 section of tlie near-arrested bicontinuous structure in the run 
of Fig. 2, showing velocity streamhnes for the two component fluids under cross-flow forcing 
with velocity of order O.Ola /t] 12j. The group of beige-colored steamlines starting at the 
bottom left show flow in one fluid (the yellow side of the interface) moving rightward; the 
second group (blue) starting on the upper right show flow in the second fluid moving to 
the left. There is no discernible motion of the interfacial structure at this flow rate. 
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Appendix: Supporting Material 

Methods 

In the simulations reported here, we choose for simphcity a perfectly symmetric pair 
of fluids with equal density p and viscosity rj. The phase diagram that controls their 
demixing is also symmetric, being described by the free energy functional (SI) 



where the order parameter Lp describes the fluid composition, and the choice of the pa- 
rameters A, B, and k controls the fluid-fluid interfacial tension a and thickness C, (SI). 

The two solid-fluid interfacial tensions are exactly equal, and the interfacial thermody- 
namics implemented as in (S2). We choose a deep quench in which the fluid-fluid interface 
should be sharp on the scale of a colloid (see below) and the thermal excitations of the 
interface (thermal capillary waves) are negligible. Thermal noise is however fully included 
in the description of fluid momentum — this imparts Brownian motion to the colloids. 
Physically relevant control parameters are then the viscosity rj, the particle radius a, the 
fluid-fluid interfacial tension a and the thermal energy ksT. (The fluid density p can be 
scaled out.) 

For our binary fluid system we use the lattice Boltzmann method, incorporating spher- 
ical solid particles (S2). Our code currently runs under OpenMP on shared memory 
parallel machines, while an MPI (Message Passing Interface) version is currently under 
development for use on larger, distributed memory machines. For the binary fluid alone, 
spinodal decomposition has been extensively studied with this code (SI) and our method- 
ology is well validated in that context. For colloids we use a standard 'bounce-back on 
links' method (S3), modified to allow for the binary solvent (S2,S4)- Brownian motion is 
achieved by a fluctuating fluid method distinct from, but closely related to, that of Ladd 




(1) 



17 



(55.56) . 

The characteristic length and time scales associated with the physics of coarsening are 

(51.57) Lq = rf/{pa) and = rf/{pa'^), which for the physical parameters chosen in 
the main text are Lq ^ 14 nm and to = 0.22 ns. Computing the same quantities in the 
'lattice units' of LB (SI) allows length and time scales to be matched to experiment, in 
principle. However, not all the dimensionless control parameters can fully be matched in 
practice. For example the Reynolds number. Re = {dL/dt)pa/ri, which characterises the 
relative importance of fluid inertia to viscosity, can be made small compared to unity, but 
not as small as the physical value (S8). We fully match e/kBT and a/LQ as described in 
the main text. The identification of our longest runs as 300 ns in duration then follows 
from the definition of to- To represent fluid-fluid-solid wetting behavior fauthfully we also 
require 1 ^ ^ ^ a in lattice units; but this is only marginally achievable (see below for 
values) . 

We performed two simulations using a lattice of 128x128x128 sites with periodic 
boundary conditions (SI). The main production runs took around 1 week on a 32- 
processor IBM p690+ system and some 2-3 weeks on a 48-processor Sun E15K system. 
These were accompanined by further runs at the same scale, and many smaller 64^ and 
32^ runs, to check that physical trends were as expected. If the system is too small, fluid 
motion will be artificially arrested once L approaches the box size (with the interface then 
attaining a state of zero mean curvature in three dimensions) even in the absence of a 
monolayer of particles ( SI ). This state was often reached with our smaller system sizes. 
However, Fig.3, which shows the domain length scale as a function of time (computed 
following (SI)), confirms that in our 128'^ runs, L remains significantly less than the 
system size: the drastic slowdown of coarsening is not a finite-size effect. 

Our free energy parameters were A = —0.002, B = 0.002, and k = 0.0014 giving an 
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interfacial thickness of ^ = 1.14, tension a = 0.0016, fluid density p = 1, and viscosity 
= 0.1 (all in lattice units). The fluid was initialised to be well mixed and at rest. A 
small amplitude random noise was added to the if field to induce spinodal decomposition. 
At the same time, colloids were positioned at rest randomly throughout the system. 
Thermal fluctuations appropriate to a temperature of 300 K were included (S5). The first 
simulation (Ref.j2ni) is a monodisperse suspension with 8229 particles of radius a = 2.3 
lattice units (corresponding to 5.4 nm in physical units) providing a solid volume fraction 
of 20%. The bidisperse simulation (Refs.|2T| l22]) has 4114 colloids of radius a = 2.3 and 
2407 larger particles of radius a = 2.74 lattice units. The overall solid volume fraction is 
again 20%. Both simulations were run initially for 520,000 time steps, which is 275 ns in 
physical time; the bidisperse run was then extended, to examine cross- flow (see below). 

Note that longer physical time scales (of order milliseconds) would be achieved if we 
chose parameters to model micron-scale colloids rather than nanocolloids. However, it 
is currently not practicable to run for timescales very long compared to the Brownian 
relaxation time of a free colloid (of any size) while still maintaining realistic values for 
Lq and to as required for the coarsening problem ( S8 ). On the other hand, the Brownian 
timescale can be reached by artificially raising the temperature; this was done for the 
structural motifs (cylinder and ripple) as described below. 

On the lattice the colloids are discrete, block-like, objects. To take account of this, a 
calibration is performed to compute an appropriate hydrodynamic radius ah (S3). This 
is the radius of the sphere which exhibits the same mean Stokes drag factor Qirrjah as 
the discrete colloid on the lattice. For the viscosity used here (77 = 0.1), the actual 
and hydrodynamic radius for the smaller colloids are found to be the same, a = ah = 2.3, 
while for the larger particles the hydrodynamic radius is slightly larger (a = 2.74 and ah = 
2.75). Fluid-mediated interactions between the particles are well represented within LB 
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when colloids are separated on the lattice, but this breaks down when the colloid-colloid 
separation h drops below the lattice scale. This can be rectified by a standard procedure 
in which lubrication forces are added by hand (S2,S3). In our runs, the normal {h~^) 
component of the pairwise lubrication interaction is corrected at interparticle separations 
h < h* = 0.7 lattice units. A much weaker transverse component of the lubrication force 
is neglected (S2). 

The computation of the lubrication forces itself becomes a major numerical exercise, 
with bad (N^) scaling in the number of colloidal particles in simultaneous mutual 
lubrication contact (S2). Since sequestration at the fiuid-fiuid interface results in very 
large N, a workround for this is essential. We achieve it by adding an additional pairwise 
thermodynamic potential (oc which effectively prevents particles approaching closer 
than roughly 0.3 lattice unit. This results in a visible residual spacing between particles 
in the interfacial monolayer (see figures). Such short range repulsions are quite common 
physically (S9) and do not seriously compromise the realism of our simulations. 

In contrast to this treatment of the lubrication forces, no equivalent corrections to the 
interparticle forces take place in the thermodynamic sector. Thus, for a dense particle 
layer, there may be relatively few fiuid-fiuid nodes left in the interstices. This under- 
estimates interfacial energies, and could account for continued slow coarsening and de- 
tachment after particles become densely packed on the interface (Fig. 3). In particular, a 
narrow neck, only one or two particles across, could become internally 'dry' with no fiuid 
nodes containing the enclosed solvent. This problem is alleviated for our higher resolution 
runs for structural motifs (Figs. SI, S2); these have somewhat larger particle size and 
considerably larger length-scales for the interfacial structures (necks, bumps) themselves. 

For the cross-fiow simulation (Fig.4), a (/^-dependent body force was applied to the 
fiuid which drives the different phases in opposite directions. A pumped flow (driven by 
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pressure gradients) is not expected to differ significantly. The body force was switclied 
on after near-arrest was complete at 520 000 time steps and tlie simulation run on to 
600 000 steps to allow a near-steady flow to establish. The mean of the velocity was 
recorded over the final 20 000 steps and used to generate the streamline ribbons seen in 
Fig.4. A small number of free particles (not attached to the interface) move discernibly 
during this time; the interfacial motion itself is negligible (the image shows the particle 
and interfacial configuration at 590 000 time steps). In a separate run (not shown) we 
increased the forcing to test the resilience of the structure. A transient (~ 0.5 ns) forcing 
twenty times stronger than in Fig.4 led to significant distortion of the interfacial structure 
followed by partial elastic recovery when forcing was removed. However, it was apparent 
that this forcing, if maintained, would lead to structural meltdown. 

To check the role of Brownian motion, we switched this off midway through a run 
whose initial parameters were that of Ref . [20] • There was a reduction in visible wobbling 
of particles at the interface, but little effect on the macroscopic motion. Reducing the 
thermal noise level is equivalent to increasing the particle radius; even with no noise from 
the outset of the run, we found very similar results to those presented above. For the 
structural motif work, we also tried artificially increasing the temperature 16-fold part-way 
through the cylinder and sinewave runs. Enhanced wobbling of particles was seen, but 
no major structural relaxation occurred on the Brownian timescale of the particles. This 
suggests that the physics of arrest is largely independent of Brownian motion and hence 
of particle size. It remains possible that structural relaxation takes place on a timescale 
well beyond the Brownian time r = Girrja'^/kBT, and thus not probed in our simulations. 
However, since the interface has come virtually to rest on this timescale, and the energy 
barriers to both particle ejection and particle rearrangement are set by avra^ ^ kBT, 
it seems unlikely that Brownian motion could overcome these barriers unaided. Similar 
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checks were made for the role of attractive bonding interactions between colloids. It is 
possible that very strong attractions (as might be required to compete effectively with 
interfacial forces) could have a strong effect, but for bonding energies of up to several 
times ksT we found no discernible difference from the purely repulsive runs described 
above. 

The higher-resolution studies (Figures S1,S2) were done using the same fluid param- 
eters as before, but with somewhat larger average particle sizes (2.7 and 4.1 lattice units 
for the cylinder, 2.1 and 3.2 for the rippled surface). These particle sizes correspond to 
physical radii of 5 nm - 9 nm for our model aqueous/hydrocarbon mixture {t] = 10~^ Pa s, 
p = 10^ kg m~^ and cr = 6 x 10^^ N m~^ at 300 K). For these studies, the thermodynamic 
interaction potential used to maintain a nonzero surface-to-surface contact distance h in 
the packed film comprised a screened Coulomb interaction with Debye length A = 0.2 
lattice units, truncated with offset so as to vanish at and beyond h = 0.4 lattice units. 
The amplitude of the interaction force (effectively, the surface charge (S9)) is chosen to 
engineer an equilibrium value of h = 1.0 lattice unit for a regular triangular lattice of 
particles of the harmonic mean size as estimated by minimizing the total energy of a unit 
cell of this lattice. This maintains surface-to-surface spacing of particles in a dense layer 
at of order one lattice spacing, which, for the chosen particle sizes, ensures that the ther- 
modynamics of the fluid-fluid interface in the interstices between particles is adequately 
resolved by the discretization. 
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Figure SI: Time evolution of a cylinder coated with bidisperse, neutrally 
wetting colloids. Left frame (particles shown translucent) shows the perturbed interfacial 
configuration shortly after initiation (5 000 timesteps). Instead of growing (as occurs for 
the particle- free cylinder via the Rayleigh-Plateau instability, culminating in pinch-off), 
the perturbation decays to a smaller amplitude and then arrests: second frame, particles 
again translucent, after 100 000 timesteps; third frame, 200 000 timesteps. There is 
almost no visible evolution between this and 600 000 timesteps (rightmost frame); without 
particles, rupture occurs at t = — 55 000 timesteps. The system occupies a 32 x 32 x 96 
lattice; particle radii are 2.7 and 4.1 lattice units. The initial state was created using a 
planar packing generated by the algorithm of (S 10) on a. periodic domain with coordinates 
then transferred to the cylinder. 
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Figure S2: Time evolution of a periodically rippled interface coated with 
bidisperse, neutrally wetting colloids. The system occupies a 96 x 48 x 48 lattice; particle 
radii are 2.1 and 3.2 lattice units. The initial state was created using a planar array of 
spheres generated by the algorithm of (SIO) on a flat periodic domain with coordinates 
then transferred to a repeating array of hemi- cylinders (upper left; after 1 000 timesteps). 
Subsequent frames (upper right, lower left, lower right) show 40 000, 80 000, and 150 000 
timesteps. Almost no macroscopic motion of the jammed monolayer occurs after t ~ 80 
000 timesteps. The evolution of a neighboring uncoated fluid surface is also shown; this 
collapses to flatness roughly as e"*/^ with r ^ 5 000 timesteps. 
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